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FOREWORD 


This  report  is  one  of  a  series  of  reports  produced  under  Office 
of  Naval  Research  Contract  No,  Nonr-4 186(00),  concerned  with  "Hydro¬ 
dynamic  Effects  Influencing  Aircraft  Carrier  Landing  Operations".  Two 
reports  have  previously  been  issued  on  the  subject  of  experimental  model 
scale  measurements  of  die  flow  disturbances  in  the  air  wake  of  an  aircraft 
carrier.  The  present  report,  dealing  with  ship  motions  at  sea,  will  be 
followed  by  a  fourth  report  of  this  series  which  is  concerned  with  time 
history  prediction  of  ship  motions. 
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ABSTRACT 


A  mathematical  model  is  developed  for  representing  the  heave,  pitch 
and  roll  motions  of  an  aircraft  carrier  at  sea.  The  data  is  in  the  form  of 
transfer  functions  relative  to  the  wavee^  which  are  determined  for  a  range  of 
forward  speeds  and  headings  considered  appropriate  to  carrier  operations 
during  the  landing  phase.  From  this  information,  spectral  responses  reprr:.senting 
the  statistical  characteristics  of  motions  in  certain  rpecific  sea  states  are  then 
developed.  The  resxilts  are  presented  in  a  form  that  can  be  applied  to  computer 
simulation  studies  for  various  types  of  wave  disturbances,  such  as  storm  condi¬ 
tions,  swells,  and  combinations  of  such  wave  systems. 

The  motions  of  heave  and  pitch  are  emphasized  in  this  study,  since  roll 
motions  are  sufficiently  small  so  that  they  do  not  significantly  influence  the  air¬ 
craft  landing  operation.  Certain  limited  conclusions  as  to  motion  severity  and 
methods  for  reducing  ship  motions  are  outlined,  together  with  recnnuuendations 
for  specific  applications  and  extensions  of  the  theoretical  results. 
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NOMENCLATURE 


Symbol 

A^tt:  ) 
e 


A 

z 


“ij 

B* 

iBM| 

CB 

CG 

c 

IGmI 

g 

h 

It 

I 

J'"e> 


Dfeflcrlption 

Spectral  energy  distribution  of  waves  referred  to 

moving  axis  system 

Component  of  added  mass  tensor 

Ratio  of  wave  amplitude  to  heave  amplitude  for  vertical 

motion-induced  wave 

Amplitude  of  surface  wave  elevation 

Matrix  element  for  vertical  plane  motions 

Local  beam  of  ship  section 

Vertical  distance  between  center  of  buoyancy  and  metacenter 
Center  of  buoyancy 
Center  of  gravity 

Propagation  speed  of  svirface  wave 

Total  energy  (integrated  eu:ea)  \mder  spectral  c\xrve 

Vertical  distance  between  CG  and  metacenter 

Acceleration  due  to  gravity 

Average  half -draft 

Total  roll  moment  of  inertia  (about  x-axis) 

Pitch  moment  of  inertia  (about  y-axis) 

Jacobian  of  freqxiency  transformation 

V 
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Symbol 

K 

KBi 


IKGi 

L 

M 


m 

o 

N' 

zz 

I  OB  I 
lOG  I 
;OMl 
p=  0 

S 


V 

V 

V 

o 

V 

w 


Deacrlption 

Roll  moment,  positive  about  the  x>axls 
Vertical  distance  from  keel  to  CB 
Vertical  distance  from  keel  to  CG 
Length  of  ship 

Pitch  moment,  positive  about  the  y-axis 
Mass  of  ship 

Local  heave  damping  coefficient 
Roll  damping  coefficient 
Vertical  distance  from  CB  to  waterline 
Vertical  distance  from  waterline  to  CG 
Vertical  distance  from  waterline  to  metacenter 
Roll  angular  velocity 

Spectral  energy  of  waves  referred  to  a  fixed  position 
Cross-sectional  area  of  ship 
Laplace  transform  op>erator 

Response  operator  (amplitude  and  phase  relative  to  wave) 
Transfer  functions  for  heave,  pitch,  and  roll,  respectively, 
for  response  to  swell 
Forward  speed 

Lateral  body  velocity  (along  y-axis) 

Lateral  wave  orbital  velocity 
Wind  speed 
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Symbol 

W 

w 

o 

X,  y,  Z 


V  ^ 

Y,Z 


7 

V 

9 

X 

P 


Description 

Ship  displacement 

Vertical  wave  orbital  velocity 

Orthogonal  axis  system;  x-aocis  horizontal  positive  toward 

the  bow;  y-axis  horizontal  positive  to  port;  z-axis 

vertical  positive  upward.  Also  used  to  represent 

translational  ship  motions  of  sway  (y),  and  heave  (z) 

relative  to  a  state  of  rest 

Bow  and  stem  x-coordinates,  respectively 

Forces  acting  positively  along  the  y  and  z-axes»  respectively 

Local  sectional  vertical  center  of  buoyancy  position 

Heading  angle  of  waves  relative  to  ship  heading 

Heading  angle  of  em  individual  wave  relative  the  wind 

Displaced  volume 

Damping  factor  for  swell  spectrum 
Elevation  of  free  siu^face 

Pitch  angle,  positive  for  rotation  about  y-axis 

Wavelength 

Fluid  mass  density 

Energy  spectrum  of  particular  motion 

Roll  angle,  positive  about  the  x-axis 

Phase  angles  of  heave  and  pitch,  respectively,  with  regard 

to  wave  reference 
vii 
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Symbol 

Description 

CO 

Circular  frequency  (rad.  /  sec) 

CO 

e 

Circular  frequency  of  encounter 

CO 

s 

Center  frequency  of  swell 

Subscripts  (unless  otherwise  defined) 

bm 

E>ue  to  body  motions  (inertial  effects) 

d 

Due  to  damping 

h 

Due  to  hydrostatic  effects 

o 

Relative  to  waterline  reference  position 

w 

Due  to  waves 

via 
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THEORETICAL  STUDY  OF  THE  MOTIONS  OF  AN  AIRCRAFT  CARRIER 

AT  SEA 

INTRODUCTION 

The  modern  aircraft  carrier  essentially  functions  as  a  floating,  mobile 
airfield  at  sea.  However,  this  mohility,  which  is  of  prime  importance  for  mili¬ 
tary  tactics,  has  certain  drawbacks  of  an  operational  nature  during  the  aircraft 
landing  phase  due  to  specific  inherent  technical  factors.  While  at  sea,  the  CEu  rier 
must  achieve  a  proper  heading  (relative  to  the  prevailing  wind)  and  also  a  specific 
forward  speed,  so  that  a  certain  minimum  wind -over- deck  air  velocity  is  main¬ 
tained  along  th'^'  anted  flight  deck.  There  is  then  a  continuous  translatory 
motion  of  the  landing  deck  as  the  aircraft  is  approaching  for  a  landing.  In  addition, 
the  carrier  operates  in  the  open  sea,  where  the  effects  of  waves  result  in  con¬ 
tinuous  irregular  oscillations  of  the  ship  in  all  six  degrees  of  freedom.  The 
varying  deck  displacements  due  to  these  ship  motions  are  significant  disturbances 
that  affect  the  pilot's  capability  in  carrying  out  a  safe  landing. 

In  order  to  aid  the  pilot  in  the  landing  operation,  certain  signal  systems 
are  used  to  provide  a  terminal  guidance  reference.  A  particxilar  type  of  system 
presently  used  on  U.  S.  Navy  carriers  is  the  Fresnel  Lens  Optical  Landing 
System  (  FLOLS),  which  makes  use  of  an  optical  image  as  an  indicator  (to  the 
pilot)  of  prop)er  glide  path  (a  description  of  the  FLOLS  and  its  operating 
characteristics  may  be  obtained,  for  example,  from  References  [  ij  and  [2j), 
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Since  the  actual  lens  system  is  located  on  the  ship,  and  thereby  partakes  of  tlie 
ship  motions,  the  beam  and  the  image  of  the  FLOLS  are  affected  by  these  motions, 
and  hence  the  lens  must  be  stabilized  against  these  motions. 

On  the  basis  of  the  important  Influence  of  ship  motions  on  the  landing  pro¬ 
cess,  as  outlined  above,  a  theoretical  study  of  the  motion  characteristics  of  air¬ 
craft  carriers  in  waves  was  carried  out  and  is  reported  herein.  The  information 
obtained  in  this  program  is  in  a  form  suitable  for  use  in  simulation  studies,  as 
a  representation  of  an  important  environmental  disturbance,  which  was  the  main 
intent  of  the  study.  In  addition,  knowledge  of  the  ship  motion  characteristics  is 
necessary  for  determining  the  required  lens  stabilization  control  rule;  it  supplies 
a  base  reference  for  determining  the  effects  of  possible  ship  motion  stabilization 
methods,  such  as  use  of  amti-pitch  fins,  nnti-roll  stabilizers,  etc. ;  and  it  also 
provides  the  basic  tools  necessary  for  development  of  a  motion  (i.  e.  time 
history)  prediction  technique. 

This  work  was  carried  out  at  (X!£ANICS,  Inc.  under  the  sponsorship  of 
Air  Programs,  Office  of  Naval  Research,  Department  of  the  Navy,  under  Contract 
number  Nonr-4 186(00). 
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ASSUMPTIONS  AND  TECHNIQUES  USED  IN  ANALYSIS 

The  equations  of  motion  used  in  this  study  will  be  formulated  according 
to  linear  theory,  thereby  allowing  a  separation  between  the  longitudinal  modes  of 
motion  (in  the  vertical  plane)  and  the  lateral  (horizontal  plane)  motions.  As  a 
further  simplification,  only  the  motions  of  the  carrier  that  will  affect  the  air¬ 
craft  in  its  longitudinal  motions  will  be  considered.  Thus  heave,  pitch  and  roll 
motions  of  the  carrier  will  be  studied  (roll  motion  is  included  since  it  influences 
the  vertical  displacement  of  the  landing  deck  and  also  the  FLOLS  lens  unit). 
Heave  and  pitch  motions  will  be  coupled,  and  surge  effects  will  be  neglected  due 
to  their  expected  small  influenoefor  the  present  problem.  The  roll  motion  will 
be  assumed  to  be  coupled  with  sway  motion,  and  effects  of  yaw  will  be  neglected. 
This  is  based  on  an  assumed  small  yaw  influence,  as  well  as  the  fact  that  yawing 
is  appreciably  influenced  by  the  helmsman's  inpmt  to  the  ship  rudder,  which  is 
an  extraneous  influence  not  directly  ascribed  to  an  environmental  (i.  e.  due  to 
wave  forces)  input.  Only  uncontrolled  ship  motions,  in  direct  response  to  wove 
effects,  are  considered  in  this  study. 

The  equations  of  motion  in  regular  waves  are  formulated  according  to 
linear  theory  by  the  balance  of  inertial,  damping,  restoring,  exciting,  and 
coupling  forces  and  moments.  Both  hydrodynamic  and  hydrostatic  effects  due  to 
the  body-fluid  interaction  are  included  in  the  analysis.  The  hydrodynamic  forces 
and  moments  such  as  damping,  exciting  effects  due  to  waves,  etc.  will  be 
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determined  by  application  and/or  extension  of  the  methods  described  in [3]  and 
[4j,  which  are  representative  of  a  technique  known  as  strip- theory  (also  referred 
to  as  a  form  of  slender-body  theory).  This  method  makes  the  assumption  that, 
for  an  elongated  body  with  transverse  dimensions  small  compared  to  its  length, 
the  flow  at  any  cross-section  is  independent  of  the  flow  at  any  other  section,  and 
hence  the  flow  problem  is  reduced  to  a  two-dimensional  problem  in  the  trans¬ 
verse  plane.  The  forces  at  each  section  are  found  by  this  method,  and  the  total 
force  is  found  by  integrating  over  the  body  length. 

Solutions  will  be  obtained  for  regular  sinusoidal  seas  with  different  wave¬ 
lengths  and  headings  relative  to  the  ship,  for  different  ship  speeds.  The  actual 
conditions  considered  will  be  forward  speeds  of  10  knots,  20  knots,  and  30  knots; 
wavelengths  extending  from  half  the  ship  length  up  to  2.  5  times  the  ship  length; 
auid  headings  ranging  from  pure  head  seas  to  seas  45^  off  the  bow.  These 
conditions  are  considered  sufficient  to  cover  most  of  the  practical  operational 
situations  associated  with  carrier  landings.  Response  amplitude  operators 
(i.  e.  transfer  function  amplitudes  relative  to  the  wave  amplitude)  are  found 
from  these  solutions,  together  with  the  phases  of  the  motions  relative  to  the 
regular  wave  system. 

The  results  for  regular  sinusoidal  waves  can  be  generalized  to  yield 
information  on  the  motions  in  an  irregular  (random)  seaway  in  terms  of  power 
spectra,  which  represent  the  relative  amplitudes  as  a  function  of  frequency. 
These  motion  power  spectra  can  be  operated  upon  to  obtain  information  on 
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average  values,  probabilities  of  relatively  high  aonplitudes  of  oscillation,  etc. 
in  the  ship  motion  time  histories,  according  to  the  methods  of  [  5j.  The  on¬ 
coming  irregular  sea  conditions  are  represented  by  a  surface  elevation  (wave) 
spectrum  as  the  input  data  for  this  process.  Three  different  fully  developed 
sea  states  were  considered  in  this  ctudy,  viz.  Sea  States  4,  5  and  6,  v/here  the 
spectral  data  for  the  different  sea  states  were  determined  from!  6j.  Similarly, 
spectral  response  data  were  obtained  for  various  swell  conditions  (i.  e.  irregular 
waves  of  varying  amplitude,  but  with  almost  constant  period)  by  developing 
narrow  band  filters  for  wave  representations,  in  conjunction  with  the  regular 
wave  response  operators.  The  particular  techniques,  results,  etc.  developed 
for  this  purpose  will  be  discussed  in  detadl  in  a  later  section  of  this  report. 

The  particiUar  aircraft  carrier  class  studied  in  detail  in  this  program 
was  the  U.  S.  S.  Forrestal  (CVA-59)  class,  which  is  representative  of  a  number 
of  large  attack  carriers.  Some  initial  consideration  was  given  to  the  evaluation 
of  motion  characteristics  of  the  U.  S.  S.  Midway  class  as  well,  but  up-to-date 
information  on  the  present  body  lines,  configurations,  etc.  of  that  ship  was  not 
readily  available  [  7  J.  Therefore  only  some  indications  of  how  to  make  use  of 
the  Forrestal-class  data  to  determine  approximate  motion  characteristics  of 
Midway-class  carriers  will  be  given  in  a  later  section  of  this  report.  A  summary 
of  the  numerical  values  of  the  important  parameters  characterizing  the  U.  S.  S. 
Forrestal,  for  which  computations  in  the  present  study  were  carried  out,  is 
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presented  In  Table  1.  Body  lines,  section  characteristics,  etc.  were  determined 
from  (8j. 
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TABLE  1 

NUMERICAL  VALUES  OF  U.  S.  S.  FORRESTAL  SHIP  CHARACTERISTICS 


Length(LBP).  ft.  990 

Beani,  ft.  129.  33 

Draft ,  ft.  35.  8 

Displacement,  tons  78,  037 

Metacentric  height,  GM|,  ft.  12.5 

Vertical  distance  from  keel  to  CG,  j  KCj ,  ft.  44.  6 

Vertical  distance  from  keel  to  CB,  j  KB|,  ft.  19.  65 

Pitch  moment  of  inertia,  without  added  fluid  inertia^,  slug-ft?  335.  97  x  10^ 

2  2  9 

Total  roll  moment  of  inertia  (including  added  fluid  inertia)  $  8lug*ft  .  21.  97  x  10 

3 

Heave  period  ,  sec.  8.  69 

3 

Pitch  period  ,  sec.  9,  10 

3 

Roll  period  ,  sec.  19.  82 


1.  Assuming  longitudinal  gyradius  ■  0.  25  L. 

2.  Obtained  from  knowledge  of  natural  roll  period,  displacement,  and  meta¬ 
centric  height. 

3.  Obtained  from  David  Taylor  Model  Basin  letter  [  9  J  and  subsequent  discussions 
with  DTMB  personnel.  All  conditions  outlined  in  Table  1  above  were  selected 
toi  conform  to  conditions  of  model  tests  carried  oi^t  at  DTMB. 


^CEANICSr. 


3. 


EQUATIONS  OF  MOTION 

The  equations  of  motion  are  formulated  relative  to  an  axis  system  whose 
origin  is  located  at  the  CG  of  the  ship.  A  right-handed  caurtesian  coordinate 
system  is  selected  with  the  axes  fixed  in  the  body,  with  the  x-axis  positive  toward 
the  bow  (in  the  direction  of  forward  motion),  the  y-axis  positive  to  port,  and  the 
z-axis  positive  upward.  These  axes  are  defined  to  have  a  fixed  orientation, 
i.  e.  they  do  not  rotate  with  the  body,  but  they  can  translate  with  the  body.  The 
ship  angular  motions  can  be  considered  to  be  small  oscillations  about  a  mean 
position  given  by  the  axes.  The  dynamic  variables  in  this  study  are  the  dis¬ 
placements  y  and  z  along  the  respective  axes,  and  the  angulai'  displacements  0 
(roll)  and  9  (pitch)  which  are  defined  as  positive  in  the  direction  of  positive 
rotation  about  the  x  and  y-axes,  respectively.  (It  should  be  noted  that  the 
present  axis  system  differs  from  that  used  by  Korvin-Kroukovsky  and  Jacobs[3j 
in  regard  to  the  positive  direction  of  pitch  angle  9.  Results  and  equations  in  the 
present  study  can  be  related  to  those  in  [  3  J  by  changing  the  sign  of  9,  positive 
pitch  moment,  etc.  in  order  to  obtain  data  for  pitch  angle  that  is  considered 
positive  bow -up,  the  usual  naval  reference  condition). 

The  hydrodynanic  forces  and  moments  are  composed  of  terms  of  inertial 
nature  due  to  body  dynamic  motions;  dissipative  terms  due  to  damping  action; 
and  exciting  effects  due  to  the  oncoming  waves.  The  effect  of  the  free  surface  is 
accounted  for  in  the  inertial  and  wave  forces  by  frequency  dependent  factors  that 
modify  the  added  masses,  and  all  couplings  of  inertial  and  dissipative  nature  are 
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y, 

included  in  the  analysis.  For  the  dynamic  body  motions,  the  slender  body  theory 
result  in  terms  of  fluid  momentum  [4 J  will  yield  a  local  sectional  vertical  force 
given  by 


dZ^  dA* 

-■  V— -  »  -  A'  -(z  -  x0  +  2V0)  +  V(i  -  x0  +  V0)  ■  ~  ■  (1 

ebe  33  dx 

where  ^  vertical  force  due  to  dynamic  body  motions;  is  the 

sectional  vertical  added  mass;  z  is  heave  displacement;  0  is  the  pitch  angle;  and 

V  is  the  forward  speed.  Integration  of  this  result  over  the  ship  hull,  with  the 

assumption  that  stern  extremities, 

results  in 

\  \  *b 

^m'-  /  /  *’33=^*“*-''/  A-jjdx-e  (2) 


The  data  for  determining  the  value  of  was  obtained  from  the  work  of  Grim  [  lOJ, 


with  A '22  defined  by 


A>  -  c 
33  8 


where  B*  is  the  local  section  beam  and  the  coefficient  C  is  represented  graphically 
as  a  fimction  of  a  dimensionless  frequency  parameter  and  two  sectional  geometric 
parameters  (beam-draft  ratio  and  section  coefficient). 

The  damping  force  in  heave  arises  due  to  energy  dissipation  associated 
with  wave  generation  by  the  ship  motions  on  the  free  surface.  The  ratio  of  the 
amplitude  of  the  heave-generated  two-dimensional  waves  to  the  heave  motion 


amplitude  of  a  section  is  denoted  as  A  ,  and  by  equating  the  energy  balance 

z 
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The  wave  exciting  force  is  obtained  by  a  combination  of  terms  due  to 
inertial  effects  (via  slender-body  theory)  and  the  prtKlominant  hydrostatic  effect 
diat  arises  due  to  periodic  buoyancy  alterations  as  the  (sinusoidal)  waves  progress 
past  the  sh4>  hull.  The  hydrostatic  restoring  force  due  to  static  vertical  dis¬ 
placement  (on  a  calm  free  surface)  of  a  ship  section  is 

-  -  pgBK*-x9)  ,  (7) 
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leading  to  a  total  force  given  by 

*b  *b 

-  pg  J  B*dx  *  z  +  pg  f  B*  xdx  •  0  (8) 

X  X 

8  8 

In  a  similar  manner*  the  buoyancy  term  in  the  wave  force  is  given*  for  a  local 
ship  section*  as 


dZ 


(1) 


-JP  •  pgB*n(x,t)  (9) 

where  ii(x*t)  is  the  wave  elevation.  The  waves  are  assumed  to  be  propagating  in 
a  direction  defined  by  the  angle  p*  where  is  the  angle  between  the  x-axis  and 
the  normal  to  the  wave  crests*  which  lies  in  the  range  ■  ~  1  ^  ^  .  The 

wave  propagation  speed  c  is  taken  to  be  <  0  for  head  seas*  i.  e.  c  ■  -  |  c|*  and 
the  expression  for  n  (x*  t)  in  terms  of  the  moving  (i.  e.  translating  in  the  x-direction) 


axis  system  is 


T)(x*t)  ■  a  sin 


r2wx 

L  X 


COSP  '*’(*>  t 

e 


(10) 


where 


b)  ■  T—  (V cosp  +  Id) 

0  A 


(11) 

(12) 


and  the  wave  elevation  is  evaluated  along  the  ship  centerline*  i.  e.  y  ■  0*  in 
accordance  with  linear  theory.  The  total  hydrostatic  wave  force  is  then 


dx  ainut  t 

e 
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The  hydrodynamic  part  of  the  wave  force  U  assumed  to  be  primarily  due  to 
inertial  effects  only,  and  is  obtained  by  application  of  slender-body  theory  [4j, 

[  llj  carried  over  to  the  surface  ship  case.  The  local  force  may  be  represented 


by 


dZ 


(2) 


:  n 


D 


where  S  is  cross-sectional  area,  the  operator  ^  is 


a4) 


-L.  ±.vi- 

Dt  at  Ox 


(15) 


and  w^  is  the  vertical  orbital  velocity  (at  the  depth  of  the  CB),  given  by 


w 


2tbIcI  ~2vh/\  f2¥X  a  .  n 
— — e  cos  [— 3^  cos  p  +  w^tj 


(16) 


The  CB  depth,  denoted  as  hin  Equation  (16),  is  included  here  as  an  approximate 
allowance  for  the  exponential  decay  of  waves.  The  fir<*l  expression  for  the  hydro- 
dynamic  wave  force,  obtained  by  evaluating  Equation  (14)  and  integrating  over 


the  ship  length,  is  given  by  r  ^ 

^  (2)  _  2Trag ^ - 2irh/ X | costu^t  J  (pS+A'^^)  sin( 

XMI  \  V  W  V  * 


2vx 


cos  p]  dx 


s 


sin 


V  [/ 


V 

■*■10)^°®^  J  "  33 

'  X 

s 


cos  pi  dx 


c 


(pS  A'gg)  cos  cosp;  dx 


CO*  p) 


s 

.  V 

■"iTT’®'--’  -33 

X 

s 


Pj  O'**]!  ' 


(17) 
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and  the  total  wave  excitation  force  is  then 


Z  *  Z  +  Z 
w  w  w 


(18) 


The  pitch  moment  terms  corresponding  to  all  the  force  terms  derived 
in  the  foregoing  are  easily  obtained  by  use  of  the  relation 


dM  ^ 

dx  "  ■  *  dx  * 

followed  by  the  required  integrations  over  the  hulL  I'his  leads  to 


(19) 


“bm*/  *’33 ’“**■  *■/  A'33*^  /  ^'33*' ■* 


s 


8 


/  A-jjdx  •  8  , 


s 


(20) 


M  .  “  r  N'  xdx  •  z  ~  f  N'  xdx'O  +  vf  N'  xdx  •  0  , 
a  J  zz  J  zz  J  zz 

XX  X 

SB  S 


(21) 


2 

“  Pg  /  B*xdx  *  *  ”  pg  J  B*x  dx  •  0  , 


s 


-  Pga 


s 


/  "b  /2  ' 

I  cos  J  xB*  sin 


dx 


(22) 


+  sin  w^t  J  xB*  cos^  -^cos  P  /  dx 


(23) 


i 


^CEANIC^ 


14. 


»(pS+A' jj)  .In/^p-cos  pj  d» 


s 


+i7r  “■  P  /  *-^’33  7  * 

X 

8 

8  CO®  ^ 


f^A' 

2tric|  ^  33 

8 


+  8in(i>  t 

e 


f  x(pS+  A'jj)  cos  cos  ^  dx 


(24) 


8 


X 


VX  /•**’. .  _,^/ 

X 

8 


*2^  /  *'33*H¥“®Pj  * 


The  equations  of  motion  are  given  by 


z.  +  z .  +  z.  +  z 

bm  d  h  w 

(23) 

M.  +  M  .  +  +  M 

Dm  d  n  w 

(26) 

where  m  ■  mass  of  ship,  I  ■  moment  of  inertia  of  ship  about  the  y-axls; 
they  may  be  written  in  the  simplified  form 
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*11^  “l2*  *  “is*  *  “l4*  “l5^  *  *16®  ■  \ 


*21*  *22*  *  *23*  *  *24*  *  *25®  *  *28®  *  “w  ' 


(27) 

(28) 


where 


*b  *b  *b 

*11  •“o'"/  A’33*''*12*/  N'**‘^:*iS-«/  B*dx; 

** 

Xfa  Xj^ 

*14’-/  *’33**‘J  *15  • '^  /  ^’33*“/  N't***^: 


>  > 

*18“'^/  B*  xdx;  »2J  •  -  J"  A'  -xdx; 

*»  *»  *• 

A 

*22  *  •  /  xdx  -  V  /  A'3j  dx;  .  -  pg  /  B*  xdx  ; 


8 


(29) 
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The  wave  force  and  moment  are  separated  into  the  forms 


Z  *  Z  cos(t)  t  +  Z  sin  w  t  (30) 

w  wc  e  ws  e 


M  M  cos  w  t  M  sin  (I)  t  (31) 

w  wc  e  ws  e 


and  this  allows  simple  computational  procedures  for  determining  solutions  of  the 
ship  heave  and  pitch  response  characteristics  (i.  e.  amplitude  and  phase)  relative 
to  the  oncoming  sinusoidal  waves. 

The  roll  motion  is  found  by  solving  die  coupled  roll- sway  equations,  on 
the  assumption  that  yawing  motion  per  se  has  small  influence  on  roll.  In  these 
two  degrees  of  freedom,  the  aj^ication  of  slender-body  theory  to  a  ship  section 
can  be  shown  to  lead  to  an  expression  for  the  inertial  sway  force  given  by 


dY 


bm 


dx 


"  ■  Dt  **'22''*  ■  Dt  **'42’’* 


(32) 


with  V  the  sway  velocity  and  p  ■  #  the  roll  angular  velocity.  Integration  over 
the  hull,  with  the  assumption  that  the  generalized  inertia  terms  ^'22  ^'4 

vanish  at  the  ends,  leads  to 


■^22^  ‘  ^2^ 


(33) 


where 


*b 

/ 


X 

s 


*’22'**  '  *42  ■ 


*b 

/  *’42 

X 

8 


dx 


(34) 
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In  the  development  of  the  sway  motion  equation,  since  there  is  no  restoring 
force  and  the  ma.jor  effect  is  inertial,  damping  in  sway  and  any  frequency 
dependence  in  the  inertial  terms  ^'22  ^'42  ^  neglected  (the  zero  fre¬ 

quency  limits  are  used).  The  lateral  wave  exciting  force  (neglecting  the  effect 
of  the  factor  e  S  is  given  by 


dY  Dv  ^ 

-sr  ■  ^"-dT 


(35) 


where 


2»a  cl  i  f  2»x 

V  - ^ -  Sinp  sin  i— r — COS3*rwi/ 

o  X  ^  \  X  •/ 


(36) 


which  leads  to 


' w  “  ’  s*“  P  /  |pS  +  A'22  (1  ■  ^  *  "e*}  *** 

*s  ^  (37) 


The  sway  equation  of  motion  is  then  given  by 

‘“o  ^  *22>  *  \2*  ■  \  • 


(38) 


and  it  remains  to  derive  the  roll  equation  of  motion  in  order  to  solve  the  coupled 


system. 


The  roll  equation  is  derived  on  the  basis  of  certain  simidifying  assumptions. 


The  inertial  roll  moment  due  to  body  dynamic  motions  is  expressed  as 


*S»n  ■  -  ^^42’  -  ^44'' 


(39) 
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where  .  is  the  total  added  roll  moment  of  inertia  due  to  the  fluid.  The  hydro* 

44 

static  restoring  force  in  roll  is  a  classical  result  of  naval  architecture,  and  is 
expressed  here  as 

Kjj  «  -  pg  vIGM|  0  *  -  W  iGM  I  0  (40) 

where  v  ts  the  displaced  volume,  (GM4  is  the  metacentric  hei£0it  (distance  be¬ 
tween  the  CG  and  the  metacenter  M),  and  W  *  pgp  is  the  ship  displacement.  The 
roll  damping  moment  is  represented  in  the  form 

Kd'-V  ' 

where  the  value  of  the  coefficient  is  a  constant  determined  from  available 
experimental  data  on  the  damping  ratio  obtained  in  a  roll-decay  model  experiment 
[  9j.  It  is  assumed  that  the  major  izifluence  on  roll  is  due  to  the  pure  roll  damping 
term,  and  other  motion  damping  terms  are  neglected.  Furthermore,  since  roll 
is  a  li^^tly  damped  motion,  the  value  of  damping  is  only  important  near  resonance, 
and  the  constant  value  obtained  in  the  roll  decay  tests  is  valid  for  this  purpose. 

The  hydrodynamic  added  masses  ^42  ^  ihe  roll  analysis 

are  obtained  by  integrating  the  sectional  values  (A'22  ^cl  A'^2)  over  the  hull, 
where  these  sectional  values  are  obtained  from  information  in  [  12j.  Particular 
care  is  necessary,  however,  in  the  determination  of  A '^2*  which  is  a  term 
referenced  to  the  CG  position,  while  the  values  in[  12j  are  appropriate  to  a 
reference  at  the  waterline.  Denoting  the  waterline  location  as  0,  and  with  the 
CG  of  the  ship  located  above  the  waterline,  the  roll  moment  about  the  CG  (as 
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used  In  this  study)  Is 

K  *  +  lOGj  Y  .  (42) 

where  CX}  is  the  distance  from  the  waterline  to  the  ship  CG.  For  the  con¬ 
tribution  due  to  sway  acceleration  v. 


K  -  -A^2  ^ 
o 

-  -<^2  *  A22>  '' 

o 

■  -A^2^  * 


with  the  local  value  A'^^  Interpreted  as 


A-^2*A’«2  *I°°I*'22- 

o 

The  value  of  A^^«  the  added  roll  moment  of  inertia.  Is  not  found  directly  but 
is  included  with  the  body  roll  moment  of  inertia  In  the  total  effective  roll 
moment  of  Inertia*  The  value  of  that  quantity  Is  determined  from  experimental 
knowledge  of  the  roll  nat\iral  period,  metacentrlc  hel|^t,  etc.  whldi  are 
supplied  by  [  9  J . 

The  roll  moment  due  to  waves  Is  determined  by  use  of  the  relation  In 
Equation (42),  viz. 


K  -K  tOGlY 
w  w  w 


o 

where  Y  Is  given  by  Equation  (37)  and  K 
This  term  Is  given  by 


o 


la  the  roll  moment 


(45) 

given  in  [  12j. 
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dK 

w 

I 

dx 


^(A« 
Dt  '"42 


2  I'B*]  ^ 

^  2  /  J  Dt 


3..  Dv 


(46) 


o 

where  v  is  given  by  EquatUm  (36),  z.  is  the  local  section  vertical  center  of 
o  b 

buoysncy  position  which  is  a  negative  quantity  in  the  present  case),  and  where 
additional  terms  involving  the  local  body  potential  have  been  neglected  for 
simplicity  and  in  expectation  of  their  small  contribution  to  the  total  wave  in¬ 
duced  roll  moment.  The  term  on  the  right  in  Equation  (46)  will  result  in  terms 
requiring  integration  over  the  hull,  with  sinusoidal  wei^ting  functions  de¬ 
pendent  on  the  wavelength,  but  it  has  been  found  that  for  the  present  case  the 
result  will  not  be  significantly  in  error  if  the  term  in  brackets  is  replaced  by 
simple  geometric  constants  of  the  ship,  so  that  the  approximations 


^  Dv  X.  Dv 

/  dT  dT 

X  X 

s  s 


V  ,„,  ®Dv  *b  bv 

/  siy)  /  S  -nfdx  (48) 

X  X 

a  8 

are  valid  in  this  case,  where  )0B|  is  the  distance  from  the  waterline  to  the  CB 
andjBMjis  the  metacentric  radius  (distance  from  CB  to  metacenter  M).  This 
approximation  results  in 


dK 


w 

o 


dx 


_D 

Dt 


Dv 

<^'42  v„)-p|OM|S-5j?  , 
o 


(49) 
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where  lOMjis  the  distance  from  the  waterline  to  the  metacenter,  and  in  accordance 
with  Eqv^t.on  (45)  leads  to  the  roll  moment  due  to  waves  given  by 


K 


2ira 
w  X 


"b  V 

Pgsinp  1gm|  f  S  cos  co®  §  +  ^ 


B 


-  g  sin  p(  1  -  cos  p)  J  A'^2  ^  ^ 


(50) 


8 


The  roll  equation  of  motion  then 


It#  N^#  +  w|Gm|#  +  . 


(51) 


where  1.  is  the  total  effective  roll  moment  of  inertia  (including  added  fluid 

X 

inertia),  and  after  eliminating  v  by  use  of  Equation  (38)  the  equation  in  terms 

of  pure  rolling  motion  alone  is 

. 

42 


^  o  22 


^^pg  sinp  |gm|  J  S  cosC^^cos  p 


2vx 


/  A  ^  1 

*  ¥  «  P  /  h  4  A-  (  I  -  ^COiWjcort 

^  o  22  X  “  '  ' 


cos  p+w  t)dx 

A  ® 


-(1  -  7^  cos  P)  /  A'  cos(^^  cos  p  +«  t)  dx 

|C|  <  A  e 


2«x 


(52) 
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SOLUTIC»l  OF  EQUATIONS  OF  MOTION 
Th«  ooiqpiUd  hMW  and  pitch  motions  given  by  ESquations  (27)  -  (31)  are 
solved  by  rsprossnting  the  responses  in  the  form 


and 


s  ■  z.oosw  t  ‘f  z^ainw  t 
1  e  2  e 


9  ■  9,cos<i>  t  ***  9»sinw  t 
1  e  2  c 


(53) 

(54) 


The  response  relative  to  a  unit  wave  hel^  is  found  in  the  form  of  an  amplitude 
and  phase  relative  to  the  wave  at  the  center  of  gravity  (the  eqtiation  reference 


position)  in  the  following  way : 


(SS) 


The  coefficients  of  the  heave  and  pitch  equations,  given  in  Equation  (29),  are 
functionally  dependent  upm  the  frequency  of  encounters)  ,  and  hence  the  system 
of  equations  is  then  a  linear  set  of  differential  equations  with  frequency  dependent 
coefficisnts.  Solution  of  these  equations,  even  with  a  computer,  is  complicated 
due  to  the  necessity  of  evaluating  the  coefficients  for  each  particular  frequency 
considered,  independent  of  the  fact  that  the  excitation  is  also  dependent  upon 
the  frequency.  This  situatioo  does  not  allow  a  simple  interpretatioa  of  the  basic 
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dynamics  of  the  ship  system  in  these  modes  of  motion,  since  a  single  diaracter- 
irtic  equation  does  not  result.  During  the  course  of  discussions  concerning  the 
use  of  the  res\ilts  in  a  simulation  study  of  landing  operations  aboard  aircraft 
carriers  at  sea,  it  was  suggested  [  13  J  that  the  equations  be  reformulated  with  the 
use  of  constant  coefricients  h  order  to  result  in  a  single  characteristic  equation 
that  determines  the  basic  ship  dynamics  for  a  partimilar  forward  speed.  This 
suggestion  is  based  upon  ttie  fact  that  aircraft  stability  and  control  techniques 
make  use  of  constant  coefficient  equations  (e.  g.  see  [  14])  even  though  there  is 
often  frequency  dependence  exhibited  in  the  stability  derivatives  obtained  from 
model  tests.  As  a  result  this  approach  was  carried  out  by  selecting  certain 
constant  values  of  the  varioxis  coefficients  in  ffie  equations  of  motion,  with  these 
constants  selected  in  such  a  way  that  the  basic  dynamic  characteristics  of  the 
system  will  be  exhibited.  Since  information  on  the  natural  periods  and  the  damping 
ratios  for  heave  and  pitch  were  available  (see  Table  1  and  [  9 J), values  of  the 
added  mass,  heave  damping,  added  pitch  moment  of  inertia,  and  pitch  damping 
were  selected  in  order  to  exhibit  some  agreement  with  model  data.  Typical 
examples  of  values  selected  for  the  constant  coefficient  approach,  shown  relative 
to  the  actual  variation  of  the  quality  as  a  function  of  frequency,  are  shown  in 
Figures  1  and  2.  The  final  coefficients  occurring  in  the  heave  and  pitch  equations 
of  motion  are  given  in  tabular  form  below  in  Table  2  for  the  various  spoeds  of 
interest . 
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‘24 


25 


'26 


TABLE  2 

Numerical  Values  of  Coefficients  in  Heave  and  Pitch  Equations 


.  252  X  xO 


9 


.  658  X  10 


. 218  X  10 


. 348  X  10 


12 


12 


12 


V 

=  0 

lO  kts. 

20  kts. 

30  kts. 

^11 

^12 

®13 

^14 

^15 

.  121  X  10® 

.  376  X  10^ 

.615  X  10^ 

.407  X  10® 

.245  X  10® 

.  357  X  10® 

.  469  X  10® 

. 581  X  10® 

^16 

.252  X  10® 

.  316  X  10® 

. 379  X  10® 

.443  X  10® 

^21 

^22 

.  407  X  10® 

.  245  X  10® 

.  133  X  10® 

.  206  X  10® 

-.915  X  10® 

.  350  X  10 


12 


349  X  10 


12 


.  343  X  10 


12 


(The  lines  across  the  page  indicate  a  constant  value  for  all  forward  speeds,  V) 

A  comparison  of  the  full  scale  values  of  the  heave  ?nd  pitch  natural  periods  and 
damping  ratios,  obtained  at  zero  speed  from  the  constant  coefficient  equations, 
with  those  obtained  from  model  experiment  (scaled  to  full  scale  values)  is  shown 
below  in  Table  3. 
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Comparison  of  Theoretical  and  Experimental  Naturad  Periods  and  Damping  Ratios 
Motion  Natural  Period  Damping  Ratio 


Theory 

Experiment 

Theory 

Experiment 

Heave 

8.  80  sec. 

8.  69  sec. 

.22 

.26 

Pitch 

8.  65  sec. 

9.  10  sec. 

.23 

.20 

In  view  of  the  fairly  good  agreement  shown  for  the  conditions  near  resonance, 
it  is  then  necessary  to  carry  out  more  detailed  computations  over  a  complete  fre¬ 
quency  range  in  order  to  see  if  this  approach  with  constant  coefficients  would  yield 
useful  results. 

An  example  of  results  computed  using  the  frequency  dependent  coefficients 
is  shown  in  Figures  3  and  4  (amplitude  ratios)  for  the  particular  case  of  20  knot 
forward  speed.  Similar  close  agreement  is  shown  for  the  case  of  phase  angles. 

It  can  be  seen  that  there  is  fairly  close  agreement  such  that  the  simpler  approach 
with  constant  coefficients  can  be  used  throughout  the  entire  range  of  conditions 
required  in  this  program. 

The  results  of  the  computations  of  amplitude  response  and  phase  angle, 
for  the  forward  speeds  of  10  knots,  20  knots,  auid  30  knots,  and  for  headings 
relative  to  the  bow  (head  seas)  ranging  from  0°  to  45°  (in  15°  increments)  are 
shown  in  Figures  5  through  It  can  be  seen  that  the  results  in  regard  to 
amplitude  and  i^ase,  when  plotted  as  functions  of  the  frequency  of  encounter 
cj  are  quite  close  for  the  different  headings,  in  spite  of  the  fact  that  it  is  known 
that  the  responses  for  the  same  wavelength  and  forward  speed  will  differ  for  the 
various  headings.  While  there  is  no  exact  agreement,  the  results  are  sufficiently 
close  over  most  of  the  frequency  range  where  significant  responses  occtir,  with 
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the  major  deviation  occurring  at  the  higher  frequencies  where  the  motions  are 
quite  small.  This  is  an  interesting  result,  whfch  may  possibly  be  due  to  the  fac 
that  the  ship  is  very  long  and  the  wavelengths  necessary  for  appreciable  motion 
(approximately  auid  longer)  will  result  in  approximately  equivalent  wave 
forces,  even  though  the  wavelengths  differ  for  the  different  headings  at  the  same 
forward  speed.  Checking  some  of  these  relationships  indicated  ttiat  such  a 
result  would  occur  in  a  number  of  instances,  and  it  is  possible  that  that  is  the 
reason  for  the  "collapsed" appearance  of  the  response  plots.  While  no  definite 
substantiating  experimental  data  was  obtained  in  detail,  some  preliminary 
checking  with  results  obtained  in  the  course  of  tests  at  David  Taylor  Model  Basin 
has  shown  fadr  agreemert  with  these  conclusions  for  the  cases  that  were 
checked.  However,  it  will  certainly  be  necessary  to  obtain  further  verification 
of  these  results  in  detail  with  the  complete  results  of  those  experiments  when 
they  are  published. 

The  results  obtained  in  the  present  study  are  still  useful,  since  the  in¬ 
tent  of  the  study  was  to  provide  data  for  simulation  purposes  so  that  environ¬ 
mental  influences  due  to  ship  motion  will  be  accounted  for  in  the  overall  program. 
In  the  course  of  the  checking  with  experimental  data  it  was  found  that  there  was 
fairly  good  agreement  in  regard  to  the  pitch  motion  characteristics,  with  the 
heave  motion  characte  :'istics  determined  in  this  program  underestimating 
the  experimental  values.  The  experimental  information  was  obtained  from  the 
limited  instances  of  examining  the  raw  data  obtained  during  the  David  Taylor 
Model  Basin  tests,  and  also  from  some  regular  wave  tests  carried  out 
at  the  Massachusetts  Institute  of  Technology  Ship  Model  Towing  Tank  (for  head 
seas)  in  the  course  of  another  eispect  of  this  program  [  15j.  The  fact  that  the 
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pitch  motion  is  represented  fairly  well,  and  the  heave  information  does  not  have 
such  good  agreement,  is  typical  of  certain  aspects  of  the  strip  theory  approach 
used  here  (see[  16j  and  [  17  J  for  a  discussion  of  the  limits  of  this  theory). 
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Computations  for  roll  motion  were  carried  out  using  Equation  (52),  where 
the  values  of  the  roll  damping  coefficient  was  obtained  using  the  Information 
pi  )vided  by  [  9  J  concerning  the  roll  damping  ratio.  Computations  of  the  roll 
amplitude  were  carried  out,  but  the  results  Indicated  extremely  small  roll 
amplitudes  for  the  frequency  range  considered  In  this  program.  The  only  per> 
ceptible  response  will  occur  in  beam  seas,  and  since  the  natural  period  of  roll 
is  so  long,  such  responses  will  not  be  considered  in  detail  since  they  do  not 
correspond  to  useful  operational  conditions  at  sea.  The  details  of  the  roll 
characteristics  will  be  presented  when  considering  the  response  characteristics 
In  random  seas,  which  will  be  treated  In  the  next  section. 
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STATISTICAL  CHARACTERISTICS  OF  MOTIONS  IN  RANDOM  SEAS 


To  characterize  the  motions  of  a  ship  in  a  random  sea,  where  the 
motions  themselves  have  a  random  nature,  the  function  known  as  the  energy 
spectrum  of  the  motion  must  ue  found.  This  spectrum  is  a  measure  of  the 
variation  of  the  squares  of  the  aunplltudes  of  the  various  sinusoidal  components 
of  the  motion,  presented  as  a  function  of  the  frequency  of  encounter  and  the 
wave  direction.  For  an  arbitrary  motion,  represented  by  the  i-subscript, 
the  energy  spectrum  of  that  motion  due  to  the  effects  of  irregular  random  waves 
is  represented  by 


I  I  2 

)  =  T,  (w  )  A  (u)  > 
i  e  I  iti  e  I  e 


(57) 


for  a  particular  fixed  ship  heading  in  a  unidirectional  sea,  where  A  (u  )  is  the 
wave  spectrum  and  |t^^  |  is  the  response  amplitude  operator  for  that  heading 
(amplitude  of  motion  per  unit  wave  amplitude).  Since  the  response  of  the  ship 
is  determined  as  a  function  of  the  frequency  of  encounter  oj  ,  it  is  necessary  to 
represent  the  wave  spectrum  also  as  a  function  of  so  that  the  total  area 
under  the  spectrum  (which  is  a  necessary  characteristic  in  determining 
statistical  information  on  response  characteristics)  can  be  easily  determined 
by  Integration  with  respect  to  u  .  The  wave  spectrum  used  in  this  study  is  the 

G 


Neumann  form  represented  by 


2  -fi  -227(v  w)' 

R^u,)  -  Cu)  ®  e 


(58) 


2  5 

where  C  is  an  empirical  constant  having  the  value  51.5  ft.  /sec.  ,  v  is  the 

2 

wind  speed  in  units  of  ft.  /  sec. ,  and  the  wave  spectrum  R  (u)  has  the  units 
2 

ft. -sec.  The  frequency  w  appearing  within  the  expression  of  the  wave  spectrum 


is  the  pure  wave  frequency,  related  to  the  wavelength  X  by 
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The  frequency  of  encounter  oj  can  be  shown  to  be  related  to  oj  by 


=  u  + - cos  fi 

e  g 


(60) 


and  it  is  necessary  when  representing  the  wave  spectrum  as  a  function  of  w  to 
present  it  in  the  form  given  by 


A^(cj  )  ■  R^(w(w  )  J  J(a)  ) 

e  e  e 

where  J(“g)  ts  the  Jacobian  given  by 

■  8u)  '  / - - 

®  \/l  +  — ^  cos  p 


(61) 


(62) 


(see  1 5  J  for  details).  It  is  possible  to  represent  the  wave  spectrum  for  a  non- 
imidirectional  sea,  allowing  for  angular  variation  (a  two-dimensional  spectnim) 
which  results  in  a  modification  to  the  basic  Neumann  formulation  given  by 


2p  -6  -2g  /(v  w)  2  ^  ^  a  ^  ^ 

-  Cw  e  "  w  cos  P  ,  -  ir  ^  ^  o 
IT  '^w  2  w  2 


Rlw)  = 


(63) 


0,  otherwise 


where  p  is  an  angle  measured  from  the  direction  toward  which  the  wind  is 
w 

blowing  (the  predominant  wave  direction).  In  that  case  the  motion  spectrum  that 
would  occur  for  a  {>articular  ship  heading  measured  relative  to  the  wind  direction 
can  be  obtained  by  integrating  with  respect  to  the  angle  p^,  but  for  simplicity 
in  this  particular  study  consideration  has  only  been  devoted  to  unidirectional  seas. 
An  illustration  of  the  effective  spectra,  including  the  Jacobian  term,  is  shown 
as  a  function  of  w  for  Sea  State  6  (defined  by  [  6j)  for  two  different  headings  of 
a  ship  relative  to  the  unidirectional  waves  in  Figures  17  and  18  for  the  various 
forward  speeds  considered  in  this  study. 
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The  spectral  technique  for  analyzing  rauidom  irregular  time  histories 
of  motion,  as  embodied  in  the  formulas  given  above,  is  applicable  to  linear 
systems  only,  since  in  that  case  a  unique  response  amplitude  operator  is  ob¬ 
tained.  The  spectral  techniques  used  herein  follow  from  linear  superposition  of 
the  responses  to  individual  frequency  components  contained  in  the  excitation, 
and  the  final  synthesis  of  the  effects  (in  a  mean-square  sense)  is  represented 
by  the  motion  spectrum. 

From  the  spectral  density  function  O.iw  )  for  a  particular  motion,  there 

I  9 

may  be  obtained,  in  principle,  all  the  statistical  or  probabilistic  properties 
possessed  by  the  random  process.  For  example,  the  total  area  under  the 

spectral-density  function  curve,  as  defined  above, 

00 

Ej  •  / 

o 

2 

is  equal  to  ,  i.  e.  twice  the  variance  of  the  ordinates  on  the  corresponding 
time-history  curve.  Here  the  ordinate  dispersion,  or  standard  deviation,  has 
been  denoted  by  o-^^,  which  is  the  root-mean-squBu:e  value  of  the  deviations  of 
the  ordinates  from  the  mean  or  average  ordinate  which  is  assiuned  to  be  zero 
in  this  study. 

Under  the  assumption  that  the  seaway  is  a  Gaussian  or  normal  stochastic 
process  which  is  exciting  a  linear  system  (in  this  case,  the  aircraft  carrier), 
the  responses  of  the  system  will  in  turn  represent  a  Gaussian  stochastic  process. 
Characteristics  of  the  motion  time  history  may  be  obtained  in  terms  of  the 
quantity  by  relating  the  behavior  of  the  envelope  of  the  record  (interpreted 
as  the  instantaneous  amplitude  of  the  time  history  curve)  to  this  quantity. 

Such  relations  are  based  on  assumed  narrow -band  behavior  of  the  energy 


%■ 
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spectxoun  [  5j,  and  yield  expressions  for  the  meeui  amplitiide  of  oscillation 
(half  the  distance  between  the  trough  and  crest  of  an  oscillation),  the  mean 
of  the  hipest  1/3  of  such  amplitudes  (known  as  the  significauit  amplitude), 
and  other  related  statistical  parameters  of  interest  for  a  specified  sea  condi¬ 
tion.  In  particular  the  relations  for  average  pitch  amplitude  and  significant 
pitch  amplitude  are  given  by 


0  '  0.  88^E- 

av.  0 


0  . 
sig. 


(64) 


Examples  of  spectral  densities  of  the  heave  and  pitch  motions  are 
given  in  Figures  19  throu^  28  ,  which  have  been  obtained  in  accordance  with 
the  relation  of  Equation  (57).  Results  were  only  obtained  for  Sea  States  5 
and  6  since  very  small  motions  (of  no  practical  consequence)  occurred  in  Sea 
State  4.  Examination  of  these  spectra  indicates  that  the  motion  characteristics 
in  a  given  sea  state  are  largest  for  the  lower  forward  speed  and  for  headings 
that  depart  from  pure  head  seas  (data  only  valid  up  to  45°  heading).  These 
conclusions  are  only  valid  for  unidirectional  irregular  seas,  and  the  variation 
with  heading  is  due  to  the  increase  of  the  effective  wave  spectra  (including 
the  Jacobian  terms)  for  the  large  heading  angle  p,  as  can  be  seen  from  a 
comparison  of  wave  spectra  in  Figures  17  and  18  .  Siniilarly  these  con¬ 
clusions  are  based  on  the  very  similar  response  amplitude  character  .sties  ob¬ 
tained  for  different  heading  conditions  in  this  theoretical  study.  Thus  the  , 

I 

i 

results  for  the  heave  and  pitch  motion  spectra  must  be  checked  against  ex-  ^ 

periinental  data  to  determine  their  validity.  However,  certain  conclusions 

obtained  from  the  theoretical  results  are  still  evident,  viz.  that  motions  In  * 
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unidirectional  waves  are  reduced  by  heading  directly  into  the  waves  at  the 
highest  possible  forward  speed.  This  conclusion  is  only  appropriate  for  the  sea 
states  considered  in  this  study,  and  it  may  also  be  modified  by  the  effects  of 
short-crestedness  of  the  waves,  thus  limiting  its  usefulness  in  determining  optimum 
heading  and  speed  for  motion  reduction.  Since  aircraft  carriers  must  operate  during 
the  landing  phase  with  a  certain  required  minimum  wind-over-deck  velocity,  the 
operating  conditions  of  forward  speed  and  local  wind  (which  is  also  related  to  sea 
state)  are  determined  from  considerations  other  than  those  dictated  by  ship  motion. 

A  sea  condition  is  usually  associated  with  a  local  wind  value,  and  hence  the  carrier 
must  experience  the  particular  environment  associated  with  that  wind  when 
achieving  the  required  wind-over-deck  velocity.  Similarly,  any  local  swell 
conditions  in  the  area  of  operation  will  also  contribute  to  the  resultant  ship  motion, 
and  hence  the  varying  complications  that  relate  forward  speed,  local  wind,  and 
local  sea  conditions  preclude  any  simple  means  of  relating  ship  motion  to  a 
p)articular  set  of  local  environmental  conditions  (sea  and  wind). 

The  results  of  computations  of  roll  motion,  as  mentioned  previously, 
indicated  extremely  small  roll  amplitudes,  and  the  limited  extent  of  coverage  of 
the  frequency  domain  (for  wavelengths  up  to  2.  5  times  the  ship  length)  shows  that 
complete  response  characteristics  could  not  be  obtained.  However,  some  attempt 
was  made  to  extend  these  curves  by  analytical  techniques  in  order  to  allow  a 
determination  of  roll  spectra.  The  form  chosen  for  representing  the  roll  response 
operator  Is  given  by 

-  «=  T.  =  lO"^  sin  p - - - 2  (65) 

’I  1  +  .  202s  +  9.951  s 
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in  transfer  function  form  where  s  *  lu  is  substituted  in  order  to  obtain  the 
frequency  characteristics.  Using  this  information,  values  of  roll  spectra 
rrere  obtained  for  a  number  of  conditions  with  larger  heading  angles,  up  to 
beam  seas,  in  order  to  illustrate  the  ai^arance  of  the  resulting  spectra. 
Some  of  these  results  are  exhibited  in  Figures  29  through  31 ,  where  it  can 
be  seen  that  Sea  State  6  results  in  a  much  larger  spectral  response  than  Sea 
State  5. 


In  order  to  acco\uit  for  the  effects  of  swells,  a  swell  spectrum  was 
arbitrarily  developed  in  a  form  appropriate  to  a  sharply  tuned  second  order 
system  with  small  damping,  with  t-  n  sn^ctrum  (in  terms  of  the  frequency 
of  enco\mter  w^)  represented  by 


A^(u)  )  » 

s  e  s 


(66) 


The  swell  spectr  ^  area  is  given  by 

OD 

E  »  r  A^(w  )dw  ■ 

s  s  e  e  8 

o 

and  since  the  significant  wave  height  is  given  by  *  2.83 this 

determines  the  value  of  the  gain  parameter  K  for  a  given  value  of  u>  (center 

8  8 

frequency)  and  (the  damping  factor  which  determines  the  shazinesa  of  the 
swell  spectrum).  In  order  to  represent  the  effects  of  heave  motion  and  pitch 
motion  in  swells,  it  was  assumed  that  the  heave  and  pitch  response  character* 
istics  were  relatively  broad  in  the  region  of  the  swell  frequency.  Thus  only  a 
constant  multiple  (chosen  from  the  response  amplitudes  at  the  particular 
swell  center  frequency)  was  selected  to  multiply  the  swell  transfer  function 
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to  represent  the  basic  transfer  fxir^ctlon  for  heave,  which  is  given  by 


T 

z 


K  (w  )  T 
z  s 


s 


K  )K 
z  s  a 


a 


(68) 


where  K  (u  )  is  the  average  heave  response  amplitude  in  the  region  near  a>  . 

Z  8  S 

The  pitch  transfer  function  for  swell  was  constructed  in  such  a  way  that  there 
would  be  aj^roximately  a  90*^  phase  difference  between  heave  and  pitch,  and  the 
expression  for  the  pitch  transfer  function  is  given  by 


T 
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W 
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S  CJ 


(69) 
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where  K-(cj  )  is  the  average  pitch  response  amplitude  in  the  region  about  cj 

w  8  S 

obtained  from  the  amplitude  response  diagrams.  An  additional  attenuating  term 
was  placed  in  the  transfer  function  denominator  to  insure  attenuation  at  the 
larger  frequencies.  On  the  basis  of  the  expression  for  the  roll  response  operator 
given  by  Equation  (65),  the  roll  transfer  function  for  response  to  swell  is  given 
by 


T.  =  T.  T 
9  0T1  a 


lO’^  K  sin  p 
s 


s 


(1  +  .202s  +  9.951s^Kl  +  2t  —  ) 

S  CJ  4 


(70) 
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These  forms  of  transfer  functions  in  the  standard  transfer  function  variable  s 
are  formulated  in  this  manner  for  aid  in  the  computer  simulation  and  analysis 
that  was  carried  out  by  Systems  Technology,  Inc.  [  18 J.  In  their  case,  using 
an  adjoint  technique  for  system  analysis,  tney  required  the  spectra  of  the  ship 
motions  to  be  represented  in  the  form  of  an  effective  shaping  filter  that  operates 
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on  a  white  noise  input  of  unit  amplitude,  and  the  equations  for  representing  the 
swell  effect,  as  shown  above,  are  appropriate  for  this  pu*  pose. 
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APPUCATION  TO  ANALYSIS  OF  AIRCRAFT  CARRIER  OPERATIONS 

Since  the  intent  of  this  study  was  to  develop  proper  representations  of 
aircraft  carrier  motions  at  sea,  in  order  to  aid  in  a  system  smalysis  of  the  air¬ 
craft  carrier  landing  operation,  plans  were  made  to  implement  an  actual 
representation  that  would  be  useful  in  the  computer  simulation.  A  set  of 
representative  severe  ship  motion  conditions  was  established  by  selecting 
certain  wave  systems  for  analysis  by  the  methods  described  in  the  previous 
sections  of  this  report.  The  sea  conditions  were  considered  to  be  made  up  of  a 
storm  sea  appropriate  to  Sea  State  6,  together  with  a  swell  corresponding  to  a 

center  frequency  given  by  w  =  0.4,  with  a  15  ft.  significant  height.  The  Sea 

s 

state  6  storm  wave  system  was  assumed  to  have  a  17  ft.  significant  ht  Ight,  and  It 
was  assumed  that  the  heading  of  the  waves  was  30^  off  the  bow.  Assuming 
Independence  of  the  sea  amd  the  swell,  the  total  sea  spectra  Is  equal  to  the  sum 
of  the  separate  spectra  of  Sea  State  6  and  the  swell,  which  results  In  the 
response  spectra  being  the  sum  of  the  spectral  responses  to  each  of  these 
separate  systems.  The  resulting  rms  value  of  a  particular  motion,  or  average 
value,  etc. ,  are  determined  from  the  square  root  of  the  sum  of  the  squares 
of  the  separate  response  characteristics  due  to  sea  and  swell  taken  alone. 

Since  the  swell  center  frequency  usually  is  lower  thaui  the  band  of  maximum 
energy  in  a  sea  state,  it  can  be  expected  that  heave  motions  will  be  larger  in 
the  swell  than  in  a  sea.  since  the  response  amplitudes  indicated  in  Figures  6,  8 
and  10  show  a  larger  heave  response  at  the  lower  frequency.  No  definite 
statement  can  be  made  with  regard  to  pitch  in  this  case,  but  it  can  be  expected 
that  results  from  both  inputs  will  be  commensurate  if  there  is  sufficient  overlap 
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of  the  spectral  bands  of  sea  and  swell  with  significant  energy  through  the  region 
of  pitch  response.  With  regard  to  roll^  it  is  expected  that  the  spectral  response 
will  be  larger  for  swell,  since  roll  is  a  sharply  tuned  motion  with  a  low 
natural  frequency.  These  combined  motions,  which  are  more  realistic  for 
operation  at  sea,  were  evaluated  for  a  Forrestal  class  carrier  at  a  forward 
speed  of  10  knots,  and  the  results  compared  favorably  with  those  obtained  from 
full  scale  data  of  the  USS  Essex  [  19j.  This  allowed  the  use  of  representative 
fits  of  appropriate  filter  forms  to  be  carried  out  by  Systems  Technology,  Inc. 
in  order  to  carry  out  their  system  analysis  of  the  aircraft  carrier  landing 
operation  (see  [  18j  for  details  of  the  simulation  and  their  results  on  landing 
operations). 

The  evaluation  of  any  spectral  response  for  a  particular  desired 

operational  condition,  either  sea  or  swell  or  their  combination,  can  be  carried 

out  in  the  manner  outlined  above  using  all  of  the  basic  data  developed  in  this 

study.  The  major  use  of  these  results  in  the  landing  analysis  has  been 

restricted  to  the  heave  and  pitch  data,  since  the  roll  effect  does  not  appreciably 

affect  terminal  landing  dispersions  I  18j.  As  discussed  earlier,  the  basic 

conditions  established  at  sea  for  a  landing  operation  require  a  specific  wind- 

rvcr-deck  velocity,  and  the  ship  must  necessarily  adapt  to  the  environment  that 

is  present  during  the  particular  time  of  the  aircraft  operations.  The  limits 

of  ship  motion  that  will  interfere  with  the  aircraft  carrier  landing  operation  are 

a  pitch  motion  of  the  order  ol  i  1.  5°  -  2°,  which  is  not  a  large  motion  per  se, 

o 

However,  for  a  ship  of  the  Forrestal  class  a  ±  2  pitch  angle  results  in  the 
stern  ramp  varying  its  up  and  down  motion  through  a  height  range  of  35  ft.  ,  and 
it  can  be  seen  that  such  a  motion  will  certainly  limit  the  effectiveness  of  aircraft 
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landing  operations.  Thus  the  ship  motion  limits  in  various  sea  states  should  be 
established  in  order  to  determine  the  influence  of  these  motions  on  the  landing 
operation,  which  is  the  intent  of  [  18j  and  subsequent  analyses. 

Another  application  oi  the  ship  motion  data  which  is  essential  for  proper 
analysis  of  motion  influence  in  the  landing  operation  is  for  determining  the 
effective  kinematics  of  the  FLOLS  beam  due  to  the  influence  of  the  various  ship 
motions,  and  also  the  required  stabilization  technique  for  the  lens  unit  itself. 

A  detailed  analysis  of  the  kinematics  is  presented  in  [  18j,  where  the  distances 
to  the  lens  location  are  referred  to  the  effective  "center"  of  pitch  and  the 
effective  "center"  of  roll.  These  locations  are  not  fixed  points  on  the  ship,  but 
vau*y  ver  a  range  of  positions  during  a  motion  '  'me  history.  Actually  the 
pitching  axis  is  supposedly  located  at  a  pomt  which  has  "zero"  acceleration, 
or  at  the  point  of  minimum  vertical  motion  [  20j.  Since  these  locations  are  not 
precise,  it  appears  to  be  more  useful  to  refer  to  a  better  reference  point,  and 
if  the  ship  motions  determined  from  the  present  analysis  are  used,  then  the 
appropriate  relerer.i  e  point  for  the  resultant  ship  motions  should  be  the  center 
of  gravity  of  the  ship.  While  that  location  may  be  difficult  to  select  also  be¬ 
cause  of  different  loadingn  on  a  ship  at  different  times,  it  is  certadnly  more 
realistic  and  there  is  a  greater  reliability  since  the  CG  location  can  be  relaxed 
to  the  location  of  the  center  of  buoyancy  (CB),  which  determines  the  ship  trim 
conditions. 

The  information  obtained  herein  is  only  appropriate  to  aircraft  carriers 
of  the  Forrestal  class,  but  it  is  possible  to  obtain  some  information  cn  other 
carriers  from  this  data,  within  certain  limits.  The  technique  for  extending 
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these  results  to  a  smaller  carrier  class  (stich  as  that  of  the  Midway  class)  is 
to  assime  that  the  responses  of  the  ships  will  be  the  same  for  the  same  ratios 
of  wavelength  to  ship  length,  i,  e.  .  This  is  then  converted  to  the  appropriate 
frequencies  of  encounter  for  the  other  class,  taking  account  of  forward  speed 
and  direction,  etc. ,  and  a  response  amplitude  and  phase  can  be  obtained  in  this 
manner.  However,  this  equivalence  of  similar  response  for  the  f  ame  values 
of  ~  will  not  necessarily  carry  over  to  other  ships  because  the  natural  periods 
of  the  motions  of  heave,  pitch,  and  roll  are  important  in  determining  particular 
ship  responses.  Since  these  periods  do  not  necessarily  scale  by  simple  linear 
relationships  in  terms  of  length  ratios,  there  may  be  some  differences  in  the 
results.  Hr  ;,  »ver,  this  approach  may  still  be  useful  as  a  first  approximation. 

In  co'iu„tion  to  the  use  of  the  results  obtained  in  this  study  for  simulation 
ixirposes  In  a  system  analysis,  there  are  a  number  of  odier  possible  applica¬ 
tions.  The  data  obtained  herein  is  necessary  for  determining  the  stabilisation 
etfort  required  if  devices  such  as  anti-pitching  fine  are  considered  to  be 
installed  on  carriers  for  the  purpose  of  ship  motion  reduction.  The  basic 
environment  data  is  important  also  for  the  actual  development  of  automatic 
landing  systems  for  carriers,  since  it  outlines  the  environmental  characteristics 
at  sea,  which  must  be  used  to  develop  specifications  for  system  performance. 
With  the  recent  installation  of  inertial  guidance  and  navigation  systems  on 
board  carrier-based  aircraft,  it  is  necessary  to  carry  out  proper  alignment 
of  the  inertial  systems  while  on  board  the  carrier.  The  information  on  ship 
motion  contained  herein  will  be  useful  for  providing  data  appropriate  to  that 
particular  operation.  Another  siipiificant  applicatioc  of  these  results  is  for 
determining  the  required  theoretical  '^wilding  blocks"for  development  at  methods 
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that  would  allow  prediction  of  the  motion  time  history.  It  is  intended  to  proceed 
with  this  last  approach  since  it  appears  to  provide  a  significant  pay-off  in  in¬ 
creased  performance  during  landing  operations  [  18j. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  main  result  of  this  study  is  a  mathematical  model  which  yields  a 
useful  representation  of  aircraft  carrier  motions  at  sea  for  purposes  of  computer 
simulation  studies,  which  are  necessary  for  a  complete  system  analysis  of 
aircraft  carrier  landing  operations  at  sea.  While  the  precise  numerical  results 
have  not  been  completely  verified  by  comparison  with  experiments,  a  limited 
comparison  with  unpublished  and  preliminary  data  exhibits  good  agreement  for 
pitch  motion  and  poorer  agreement  for  the  case  of  heave.  Since  the  pitch 
motion  results  in  larger  vertical  deck  excursions,  which  are  more  significant 
in  the  carrier  landing  problem,  the  data  developed  herein  provides  useful  in¬ 
formation  on  the  most  significant  aspect  of  ship  behavior  at  sea  for  this 
particular  problem. 

With  regard  to  the  detailed  results  of  the  dieoretical  computations,  the 
frequency  response  characteristics  (amplitude  and  phase)  using  the  constant 
coefficient  model  are  sufficiently  close  to  the  responses  obtained  by  use  of  the 
standard  frequency- dependent  coefficient  technique.  This  approach  allows 
computational  simplicity  and  more  insight  into  the  basic  ship  dynamic  character¬ 
istics  (roots  of  the  system  characteristic  equation  which  exhibit  the  natursil 
frequencies  and  damping  ratios  for  coupled  motions).  Another  feature  of  the 
results  is  limited  dependence  on  heading  angle  of  the  frequency  response 
characteristics  of  a  large  aircraft  carrier,  when  presented  as  a  function  of  the 

frequency  of  encounter  co^.  This  result  may  have  significance  with  regard  to  the 
influence  of  wave  directionality  properties  on  the  ship  responses.  In  view  of  these 
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various  findings,  it  is  recommended  that  these  results  be  checked  with  experi  - 
mental  data  when  such  data  becomes  completely  available.  It  appears  plausible 
that  the  present  theory  (if  it  is  fairly  representative  of  actual  experimental  values) 
can  be  used  in  conjunction  with  experimental  data  to  provide  guidance  in  deter¬ 
mining  the  response  characteristics  c.  an  aircraft  carrier  for  a  large  range  of 
operating  conditions.  This  infor maticn  should  be  applied  to  determine  the  trends 
in  the  significant  variables  (such  as  sr^ed,  heading  angle,  etc. )  that  will  yield 
minimum  ship  motions  for  other  c^a  conditions  than  those  considered  in  the 
present  study.  In  particu’.  ir,  this  information  should  be  applied  to  real  spectra 
that  have  been  measured  at  sea  ^for  more  severe  conditions)  in  order  to  determine 
operational  limits. 

The  information  available  in  the  present  report  on  ship  response 
characteristics  in  waves  should  be  used  for  developing  techniques  for  motion 
time  history  prediction.  Since  transfer  functions,  spectral  data,  etc,  are  basic 
elements  in  prediction  theory,  this  appears  to  be  a  natural  extension  of  the 
results  contained  in  this  repoii.  It  is  recommended  that  an  investigation  of  the 
various  techniques  for  time  history  prediction  of  ship  motion  be  carried  out,  with 
attendant  model  tests  in  order  to  check  the  effectiveness  of  the  various  approaches. 
The  prediction  time  requirements  should  certainly  be  related  to  the  minimum 
time  required  for  specific  aircraift  maneuvers  during  landing.  Similarly,  it  is 
important  to  develop  some  means  of  ’mplementatlon  of  the  prediction  information 
into  the  landing  system  procedure  to  obtsin  the  maximum  benefits.  It  is  ex¬ 
pected  that  such  a  program,  i(  successful  in  attaining  accurate  prediction  for  the 
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proper  time  interval,  will  yield  significant  Improvements  In  the  safety  and  per¬ 
formance  of  carrier  landing  operations. 
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Fig.  1  Added  mass  as  a  function  of  frequency 
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Fig.  3  Comparison  of  pitch  amplitudes  using  frequency 
dependent  coefficients  and  constant  coefficients 
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Fig.  4  Comparison  of  heave  amplitudes  using  frequency 
dependent  coefficients  and  constant  coefficients 
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Fig.  5  Pitch  amplitude  variation  with  encounter  frequency 

and  heading,  V  =  10  knots 
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Fig.  6  Heave  amplitude  variation  with  encounter  frequency 

and  headii  g,  V  »  10  knots 
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Fig.  7  Pitch  amplitude  variation  with  encounter  frequency 

and  heading,  V  20  knots 
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Fig.  8  Heave  amplitude  variation  with  encounter 
frequency  and  heading,  V  =  20  knots 
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Fig.  9  Pitch  amplitude  variation  with  encounter  frequency 

and  heading,  V  ■  30  knots 
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Fig.  10  Heave  amplitude  variation  with  encounter  frequency 

and  heading,  V  ■  30  knots 
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Fig.  12  Variation  of  heave  phase  angle  with  encounter  frequency  and 

heading,  V  =  10  knots 
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Fig.  13  Variation  of  negative  pitch  (-  0)  phase  angle  with  encoxinter 

frequency  and  heading,  V  =  20  knots 
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Fig.  18  Effective  wave  spectrum  for  unidirectional  Sea  State  6, 


Fig,  17  Effective  wave  spectrum  for  unidirectional  Sea  State  6, 


rad.  -sec 


V  =  10  knots 


Fig.  10  Pitch  spectrum  in  unidirectional  Sea  State  5,  V  =  10  knots 
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21  Pitch  spectrum  in  xmidirectional  Sea  State  5,  V  =  30  knots 
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Fig.  22  Pitch  spectrum  in  unidirectional  Sea  State  6 
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Fig.  23  Pitch  spectrum  in  unidirectional  Sea  State  6, 

V  =  20  knots 
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Fig.  27  Heave  spectrum  in  unidirectional  Sea  State  6,  V  =  20  knots 
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Fig.  29  Roll  spectrum  in  unidirectional  Sea  State  5, 

V  =  20  knots 
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Fig.  31  Roll  spectrum  in  unidirectional 


Sea  State  6,  V  =  10  knots 
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